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Graphene exhibits unusual electronic properties, caused by a linear band structure near the Dirac point. 
This band structure is determined by the stacking sequence in graphene multilayers. Here we present a novel 
method of microscopically controlling the band structure. This is achieved by epitaxy of graphene on 
3C-SiC(lll) and 3C-SiC(100) thin films grown on a 3D microfabricated Si(lOO) substrate (3D-GOS 
(graphene on silicon)) by anisotropic etching, which produces Si(lll) microfacets as well as major Si(lOO) 
microterraces. We show that tuning of the interface between the graphene and the 3C-SiC microfacets 
enables microscopic control of stacking and ultimately of the band structure of 3D-GOS, which is typified by 
the selective emergence of semiconducting and metallic behaviours on the (111) and (100) portions, 
respectively. The use of 3D-GOS is thus effective in microscopically unlocking various potentials of 
graphene depending on the application target, such as electronic or photonic devices. 



Carrier dynamics in monolayer graphene, a honeycomb lattice consisting of carbon atoms, are described in 
terms of quantum electrodynamics^"^, as a consequence of the linear band dispersion relation. Graphene 
therefore possesses excellent electronic properties, such as giant carrier mobilities. In addition, multilayer 
graphene has the unique feature that its electronic properties, e.g., the linearity of the band dispersion, are tunable 
by changing the stacking sequenced Graphene is a promising material for various next- generation device appli- 
cations, such as high-electron-mobility transistors (HEMTs)^'^, saturable absorbers^'^ and terahertz lasers^^. 

One of the major challenges to such device applications of graphene is that of establishing epitaxial growth of 
graphene on practical substrates using methods compatible with existing device technologies^. The chemical 
vapor deposition method (CVD) produces large-scale, high-quality graphene^ but it also requires that a 
transfer process be performed to move the graphene from a metal substrate to an insulating substrate, during 
which it is hard to completely avoid metal contamination^^. Epitaxial growth of graphene on SiC bulk crystals also 
produces large-scale, high-quality graphene on SiC substrates, which are already in use for power device applica- 
tions^^'^^. Epitaxial graphene (EG) has, however, a major drawback: the high production cost of the wafer. To 
overcome this drawback, epitaxial growth of graphene on Si (GOS) substrates by sublimating surface silicon 
atoms of 3C-SiC thin films on the Si substrates at elevated temperatures has been developed, up to the wafer 
scale^^"^^ Formation of graphene has been confirmed by Raman spectroscopy^^"^^ low-energy electron diffraction 
(LEED)^^"^° and photoelectron spectroscopy (PES)^^'^°. Despite the mediocre quality of initial GOS, recent 
improvements in GOS technology now provide a material that clearly demonstrates the linear dispersion of 
the n band near the Dirac point by angle- resolved PES^\ The potential of GOS for electronic and photonic device 
applications has been demonstrated by fabricating transistors^^'^^ terahertz photonic devices^^ and even logic 
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Figure 1 | 3D-GOS Fabrication, (a) Schematic of the 3D-GOS fabrication procedure, (b) Optical microscope image of the 3D-GOS. 



inverters^^, a fundamental component for integrated circuits. GOS is 
thus becoming a promising graphene production method, based on 
its advantages, such as the cost of the wafer and compatibility with 
existing Si-based device technologies. 

In addition to these manufacturing advantages, GOS has another 
more profound advantage: the tunability of the band structure by 
altering the Si substrate. By controlling the crystallographic orienta- 
tion of the Si substrate, we can select the surface orientation of the 
3C-SiC thin film, which eventually makes it possible to modulate the 
stacking of graphene and its interface structure with SiC. The latter 
controls the presence or absence of the interfacial (6^3 X 6^3) -recon- 
structed buffer layer, which is a precursor for graphene in EG and 
GOS. On the Si(l 1 1) substrate, the buffer layer is present, and GOS is 
Bernal-stacked^^'^^; the energy band then loses its linearity. It is, 
however, multiply- split, and a bandgap can be opened, which is 
suitable for electronic device applications. On the other hand, on 
Si(llO) and Si(lOO) substrates, the buffer layer is absent, and GOS 
is not Bernal-stacked'^''^'^^. The GOS on Si(llO) forms turbostratic 
stacks'^ while GOS on 3C-SiC(100)/Si(100) has rotational stacking 
faults with a finite angle of 15 degrees^^. The energy bands in these 
cases are anticipated to maintain their linear natures; the bands are 
not split, and the bandgap is closed. This feature is promising for 
photonic device applications, especially in the terahertz region. Such 
controllability over physical properties of graphene simply by chan- 
ging the crystallographic orientation of Si substrates is definitely one 
of the major advantages of the GOS technology. It is not easy to take a 
similar strategy by using SiC bulk crystals. If we can vary the crys- 
tallographic orientations within an identical Si wafer, however, this 
would even enrich the applicability of GOS. This is actually possible 
by use of the well-matured Si microfabrication technology. The ques- 
tions, then, are (1) whether similar graphene structures are repro- 
duced on such microfabricated, neighboring (100) and (111) 
portions as on independent wafers, and (2) whether the controllabil- 
ity of the physical properties of graphene is microscopically realized 
on (100) and (111) portions on a single Si substrate. The question (2) 
makes sense when we recall the fact that a continous few-layer gra- 
phene grows over both the (100) and the (111) portions^^'^^. 

In this work, we demonstrate the microscopic tuning of the band 
structure of epitaxial graphene on a 3D microfabricated Si substrate 
(3D-GOS) through stacking with the aid of microelectromechanical 
system (MEMS) technology, specifically anisotropic wet etching 
which produces (111) and (100) microfacets on a Si(lOO) substrate. 



Results 

3D-GOS formation. The formation of 3D-GOS is schematically 
shown in Figure 1. First, a sacrificial Si02 film (90 nm) is grown 
on a Si(lOO) substrate (p-type, 1-10 Q cm) by the dry oxidation 
method. Electron-beam lithography is performed on the sacrificial 
Si02 thin film, followed by fast atomic beam etching using SFe as 
an etchant^^'^^, to leave behind the Si02 mask on the surface. 
Using the Si02 mask, the surface is anisotropically etched in a 25% 
tetramethyl ammonium hydroxide (TMAH) aqueous solution at 
353 K for 6 min. The etching by the TMAH aqueous solution is 
used to expose Si(lll) microfacets on Si(lOO) substrates, which are 
frequently used as the standard fabrication procedure for MEMS 
devices^^. The etching depth is ~1 jam. Si(lll) microfacets then 
form on the Si(lOO) substrate, followed by the removal of the Si02 
mask with a treatment in a 5% HF aqueous solution. The 3C-SiC thin 
films (—100 nm in thickness) are grown on microfabricated Si(lOO) 
substrate by using gas-source molecular beam epitaxy with mono- 
methylsilane as a gas source at a substrate temperature of 1353 K^^ 
Due to its almost conformal growth, the 3C-SiC(lll) and 3C- 
SiC(lOO) microfacets are formed thereon. Finally, the substrate is 
annealed in vacuo at 1523 K for 30 min to grow graphene epita- 
xially on the surface of the 3C-SiC(lll) microfacet^^"^^'^^. Various 
kinds of micropatterns (Fig. lb) are formed on the sample to calibrate 
the device fabrication for future fabrication projects. 

Microscopic variation of the interface of 3D-GOS. The interface 
between graphene and the 3C-SiC microfacets is examined here 
because the graphene growth mode is determined by the presence 
or absence of the buffer layer because it works as the template for the 
Bernal stacking of graphene^^'^^'^°. The interface chemical composi- 
tion is probed nanoscopically by using three-dimensional high- 
resolution scanning photoelectron microscopy (3D nano-ESCA 
(electron spectroscopy for chemical analysis)) using a focused 
incident X-ray beam with a diameter of 70 nm and a photon 
energy of 1000 eV^°'^\ The energy resolution of the spectrometer 
was set to 300 meV^^. Figure 2a shows the photoelectron intensity 
mapping of the Cls core level (—284 eV) due to the graphene of 3D- 
GOS. The sample surface is found to be covered with graphene on 
both the (100) and (111) facets, although the image contrast 
reflecting the microfabrication pattern of 3D -GOS arises from the 
difference in relative angle of the microfacets to the detector used 
for 3D nano-ESCA. Figure 2b,c shows pinpoint Cls core-level 
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Figure 2 | Interface analysis of 3D-GOS. (a) Intensity mapping of the 3D -GOS by using CI s core-level photoelectrons (—284 eV). (b) Pinpoint CI score- 
level spectrum of graphene on the 3C-SiC(100) microterrace. (c) Pinpoint Cls core-level spectrum of graphene on the 3C-SiC(lll) microfacet. (d) X- 
TEM image of the interface between graphene and the 3C-SiC(100) microterrace. (e) X-TEM image of the interface between graphene and the 
3C-SiC(lll) microfacet. 



photoelectron spectra of graphene on the 3C-SiC(100) and 3C- 
SiC(lll) microfacets using 3D nano-ESCA. The peaks due to the 
3C-SiC thin films (-283 eV)^"'^' and graphene (-284 eV)^"'^' are 
visible in both of the spectra. The ratio of the intensity of the SiC peak 
to the intensity of the graphene peak is decreased by increasing the 
photoelectron emission angle from the surface normal (surface- 
sensitive mode), as seen in the angle-resolved spectra (see Fig. SI). 
This confirms that graphene is present on the 3C-SiC microfacets. 
The estimated layer numbers of graphene on the 3C-SiC(100) and 
3C-SiC(lll) microfacets are 2.7 and 2.4, respectively, using the 
standard equation using the intensity ratios of the graphene and 
SiC peaks (see Supplementary Note 1). The difference between the 
spectra is the presence of the shoulder in the spectrum of graphene 
on the 3C-SiC(lll) microfacet, labelled as peak B at the higher 
binding energy (—285 eV). The analysis of the angle- resolved 
spectra clarifies that peak B is ascribable to the buffer layer situated 
between the graphene and the 3C-SiC(lll) microfacet^^'^^ (see 
Supplementary Note 2). The buffer layer is thus formed selectively 
between the graphene and the 3C-SiC(lll) microfacet. A reason for 
this selective formation is that the buffer layer is commensurate with 
the 3C-SiC(lll) surface and energetically favorable^^'^^, while the 
twofold symmetry of the 3C-SiC(100) surface may inhibit the 
formation of such a buffer layer. 

The atomic structure of the interfaces between graphene and the 
3C-SiC microfacets has been investigated by using cross-sectional 
transmission electron microscopy (X-TEM), as shown in Figure 2d,e. 
Planar-shaped graphene and 3C-SiC are atomically resolved in both 
of the X-TEM images. The interlayer distance between the graphene 
layers is 0.34 ± 0.07 nm, which agrees well with the interlayer dis- 
tance of graphite as well as that of GOS using flat Si substrates^^. The 
layer numbers (2 to 3) obtained from the X-TEM images are con- 
sistent with those estimated from the pinpoint spectra. The differ- 
ence is the distance between the nearest layer and the 3C-SiC 
topmost layer. The distance (0.26 ± 0.07 nm) between the nearest 
layer (buffer layer) and the 3C-SiC(lll) microfacet is shorter than 



that between the nearest layer (graphene) and 3C-SiC(100) (0.34 ± 
0.07 nm). This shortened distance is explained by the presence of 
covalent bonds in the buffer layer that are tied to the surface of the 
3C-SiC(l 1 1) microfacet^^. No such covalent bonds exist between the 
nearest graphene layer and the surface of the 3C-SiC(100) microfa- 
cet^^. The X-TEM observation corroborates the 3D nano-ESCA 
measurement which shows that the interface varies with the crystal - 
lographic orientation of 3C-SiC microfacets grown on the microfab- 
ricated Si(lOO) substrates. 

Microscopic variation of the stacking of 3D-GOS. To clarify the 
relation of the observed interface variation caused by microfacets 
with the graphene stacking, microscopic LEED (|i-LEED) analysis 
is performed for 3D-GOS, as shown in Figure 3. The LEED images 
were taken by collecting diffracted electrons from selected areas of 1- 
jim diameter^^'^^'^^. The distortion of the |a-LEED images is due to the 
inhomogeneous surface electric field arising from the inclination of 
the (111) microfacets. This is because the diffracted electrons are 
having lower kinetic energies (—50 eV) near the surface and are 
thus susceptible to the inhomogeneity of the surface electric field^^. 

As seen from the ji-LEED images of the two bevels, a hexagonal 
LEED pattern from graphene is clearly identified on the 3C-SiC(l 1 1) 
microfacets. The spots due to the specular reflection, labelled as (00) 
in the images, are located away from the center of the |i-LEED pat- 
tern because of the inclination of the 3C-SiC(lll) microfacets. The 
|i-LEED pattern is almost the same as that of the graphene on the flat 
Si(lll) substrate and graphene on 6H-SiC(0001)^''^'''^ which indi- 
cates the Bernal stacking of graphene on the 3C-SiC(l 1 1) microfacet, 
although there can be a slight imperfection in the stacking 
sequence^^'^^. In the |i-LEED images, the spots due to the 6^3 X 
6^3 -reconstructed buffer layer are hardly visible. This can be related 
to the existence of the graphene overlayers (two to three atomic layer 
thick) on the buffer layer, which could substantially weaken the 
diffraction spots from the buffer layer^^ because of the short escape 
depth (<1 nm) of the diffracted electrons having low kinetic ener- 
gies of 50 eV'^ 
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Figure 3 | Stacking variation of 3D-GOS. |a-LEED analysis of 3D-GOS from the regions of the top (100) microterrace, bottom (100) microterrace and 
(111) bevel. The green circles indicate (1X1) spots from the Bernal- stacked graphene. The pink circles indicate the adjacent spots from rotationally- 
stacked graphene. The energy of the incident electrons is 50 eV of the |a-LEED observation. 



In contrast to epitaxial graphene on the 3C-SiC(lll) microfacet, 
the LEED spots of graphene on both the top and bottom 3C-SiC(100) 
microterraces are rotated, as observed on the graphene on 3C- 
SiC(lOO) thin films on flat Si(lOO) substrates. The rotation angle 
between the adjacent spots indicated by the yellow arrows is about 
15 degrees. One of the possible mechanisms for this rotated growth of 
graphene on the 3C-SiC(100) face is the {111} facet-induced mech- 
anism^^. In this mechanism, graphene growth on 3C-SiC(100) sub- 
strate is initiated epitaxially at the Si-terminated 3C-SiC(lll) 
microfacets, and the graphene extends over in a carpet-like manner 
toward the 3C-SiC(100) terrace. This mechanism accounts for the 
absence of the buffer layer beneath the graphene on 3C-SiC(100)^^ 
and the 15 degrees between the spots as well. The graphene's prin- 
cipal axis of < 1120 > is 15 degrees rotated from the <010> axis of 
the 3C-SiC(100) terrace in this mechanism (Fig. 3(b)). There are four 
3C-SiC{lll} microfacets surrounding the 3C-SiC(100) microter- 
race. Only two out of the four, however, are Si-terminated. The other 
two are C-terminated, which results in turbostratic stacking and a 
ring-like LEED pattern^^'^°. To explain the 24 spots using this model, 
therefore, we should consider presence of rotated domains, most 
probably of anti-phase domains. It is thus demonstrated for 3D- 
GOS that the stacking varies microscopically in accordance with 
the variation of the interface structure between graphene and the 
3C-SiC microfacet and microterrace. 

Microscopic control of the band structure of 3D-GOS. The band 
structure of 3D-GOS is anticipated to vary with the microfacet and 
the microterrace orientation owing to a high susceptibility of the 
band structure to the stacking sequence, by analogy with the 
susceptibility of the band structure of GOS using flat Si substrates 
for stacking depending on the crystallographic orientation of the Si 
substrate^^'^°. To provide proof of this assumption, the band structure 
of 3D-GOS is microscopically investigated by Raman microscopy, as 
shown in Figure 4. The fundamental vibration modes of graphene, G 
(-1600 cm"^), D (-1360 cm"^), and G' (-2700 cm"^) bands, are 
visible in the spectra on the graphene on both 3C-SiC(100) 
microterraces and 3C-SiC(lll) microfacets. The G band is the 
band originating from a single- resonant Raman process at the T 
point^^. The D band comes from double-resonant Raman 
processes, and is associated with the presence of defects^^. The 
appearance of the D band therefore indicates the existence of 
defects in graphene on the microfacets and microterraces. The G' 



band results from double-resonant processes involving strong 
electron-phonon coupling^^'^^. The G' band is the overtone of the 
D band, but is not associated with the presence of defects; rather it is 
associated with the band structure of graphene^^. The appearance of 
these modes corroborates the formation of graphene on the 
microfacets and the microterraces. 

Among these bands, the G' band is analysed in detail to clarify the 
band dispersion^^"^^ The G' band is a second-order process related to 
a phonon near the K point in graphene, activated by double-res- 
onance (DR) processes^^'^°, which are responsible for its dispersive 
nature with the excitation energy and which cause a strong depend- 
ence of the lineshape of the G' band to the band structure of gra- 
phene^^'^^'^^. As demonstrated in Figure 4, the G' band of graphene on 
the 3C-SiC(100) microterrace is not split and can be represented by a 
single, symmetrical Lorentzian. The full width at half maximum 
(FWHM) of the G' band is 70 cm~\ This value is larger than that 
of epitaxial graphene on 6H-SiC(0001) (37 cm"^)^^ because of the 
presence of the larger number of defects also inferred from the larger 
D band^\ The FWHM in this work is smaller than the previously 
reported value (80 cm"^) for graphene on a 3C-SiC(100) thin film on 
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Figure 4 | Band structure variation of 3D-GOS. Raman microscopy from 
graphene on the 3C-SiC(100) microterrace (top) and the 3C-SiC(lll) 
microfacet (bottom). The thin lines indicate the decomposed peak for the 
G' bands. The black line in the spectrum of graphene on the 3C-SiC(l 1 1) 
microfacet indicates the synthesized curve. 
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a Si(lOO) substrate'*^, which is corroborated by the larger D band in 
the previous work^^. The non-spUtting of the G' band suggests that 
the graphene on the 3C-SiC(100) microterrace has a metalUc nat- 
^j.gi7,2o,44,45 -pj^jg arises from the negUgible interlayer interaction since 
the graphene on the 3C-SiC(100) microterrace is not Bernal-stacked, 
as indicated by the |i-LEED (Fig. 3). This non-Bernal stacking is 
consistent with the absence of the buffer layer between graphene 
and the 3C-SiC(100) microterrace, which was confirmed by 3D 
nano-ESCA (Fig. 2). 

On the other hand, the G' band of the graphene on the 3C- 
SiC(lll) microfacet is asymmetrical and is much broader 
(90 cm"') than that on the 3C-SiC(100) microterrace (70 cm"'). 
These two features can only be accounted for by considering a set 
of multiple Lorentzian components, as shown in Figure 4. If the G' 
band were to be from a single Lorentzian, its FWHM should have 
been much smaller than what we observe here. This is because the 
intensity ratio Id/Ig of the D band to the G band, which is a good 
measure of the defect density that broadens the G' band, is substan- 
tially smaller in the graphene on 3C-SiC(lll) microfacet (Id/Ig — 
0.64 ± 0.04) than that on the 3C-SiC(100) microterrace (0.75 ± 
0.04). Thus, the G' band on 3C-SiC(lll) cannot be represented by 
a single Lorentzian. The broadening of the G' band of the graphene 
on the 3C-SiC(lll) microfacet is thus concluded to be due to pres- 
ence of multiple components. The presence of multiple Lorentzian 
components in the G' band suggests presence of multiple routes for 
this double-resonant Raman scattering process in the graphene on 
the 3C-SiC(lll) microfacets. Several origins can be considered, 
which include the band splitting of graphene'*^"^^, strain variation^^, 
and carrier doping variation^^ '*^. If strain is the cause, not only the G' 
band but also the G band should show a corresponding set of mul- 
tiple components'^. The observed G band, however, shows a single 
component, as shown in Figure 4. The shoulder at around 1620 cm"' 
is attributed to the so-called D' band, induced by defects^^. Therefore, 
the strain is not the cause of the multiple components of the G' band. 
Carrier doping can also shift the G' band^^'^^. In this case, again, the G 
band should also show multiple components. Moreover, the G band 
shift due to carrier doping should be larger than that of the G' band^^, 
which suggests appearance of even distinct multiple components in 
the G band. This is in contrast to the experiment. Thus, the charge 
density variation is excluded as the origin of the multiple compo- 
nents of the G' band. We therefore suggest the band splitting of 
graphene as the most likely origin of the multiple components in 
the G' band on the 3C-SiC(lll) microfacets. This assignment is 
consistent with the Bernal stacking^^'^^ revealed by |i-LEED (Fig. 3) 
as well as with the presence of the buffer layer confirmed by 3D nano- 
ESCA (Fig. 2). There is thus a good reason to expect a semiconduct- 
ing nature in the graphene on the 3C-SiC(lll) microfacets. 

Discussion 

This observed variation of the band dispersion is also corroborated 
by examining the area intensity ratio of the G' band to the G band 
because this exhibits the degree of interlayer interaction^^'^^. The area 
intensity ratio is 1.2 for graphene on the 3C-SiC(l 1 1) microfacet, and 
1.7 for graphene on the 3C-SiC(100) microterrace. The reduced area 
intensity ratio of graphene on the 3C-SiC(lll) microfacet is 
explained mainly by a stronger interlayer interaction due to Bernal 
stacking of graphene on the 3C-SiC(l 1 1) microfacet, as in the case of 
COS using flat Si(lll) substrates. In this way, the band structure of 
the 3D -COS on the microfacet and the microterrace is surely tunable 
microscopically by varying the graphene stacking. 

In conclusion, we have succeeded in microscopically controlling 
the band structure of 3D-GOS through tuning of the interfaces 
between graphene and the 3C-SiC microfacets and microterraces 
with the aid of the epitaxy on the microfabricated Si(lOO) substrate 
by anisotropic wet etching, which produces Si(lll) microfacets as 
well as Si(lOO) microterraces. The realization of 3D-GOS is the first 



step to exploiting graphene-based multifunctional integrated circuits 
(combining electronics and photonics) compatible with existing Si- 
based electronics in the next generation of devices. 

Methods 

Microscopic characterization. The spatially- resolved and angle-resolved Cls core- 
level photoelectron spectra are taken by using the 3D nano-ESCA system installed at 
the BL07LSU at SPring-8^°"^^, where the synchrotron radiation (SR) beam has a high 
energy- resolving power (£/A£ > 10*). The photon energy of the SR beam for the 
measurement is 1000 eV. Using 3D nano-ESCA enables us to obtain a high lateral 
resolution (70 nm) by focusing X-rays using a Fresnel zone plate. The energy 
resolution of the spectrometer was set to 300 meV, and the accuracy of the angle 
resolution is 0.9°. The binding energy scale was calibrated by the photoelectron peaks 
of a gold foil (Au 4/7/2). Details of the experimental setup are described in the 
previous reports^°'^\ The stacking of graphene is probed by microscopic low-energy 
electron diffraction (|i-LEED) using electron optics through a low-energy electron 
microscopy (LEEM) system installed at BL17SU at SPring-S'"'^". The |a-LEED is 
acquired by collecting diffracted electrons from selected areas of l-\im diameter. 
Raman spectra are obtained by using Raman microscopy with an excitation energy of 
2.41 eV and a lateral resolution of 1 |j,m^^. The X-TEM images are taken by using 
JEOL JEM-2010 (JEOL Ltd., Japan) with an acceleration voltage of 200,000 eV. The 
instrumental resolution of the transmission electron microscope is ±0.07 nm. 
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